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I. INTRODUCTION 



In the J-matrix approach [l,|2| to the many-channel problem the partial wave potentials 
yi a P) j n th e interaction 

? = J>)v<»0</?| (i) 

a,f3 

are given by the expansion 

iV-1 

VM(r, r>) = to £ € a (*) K,i </>%')■ (2) 

n, n'=0 

Here, a is the channel index that contains the orbital angular momentum £ a , x = r/p, the 
relative coordinate in units of the oscillator radius p = yH/ pu, where p is the reduced mass. 
In the case of collisions with neutral targets, the harmonic oscillator basis functions 



2n\ 



pT(n + £ + 



x 



(3) 



are applied. 

In [3|, |j, 15J an inverse scattering formalism within the J-matrix method has been proposed 
and developed, where potentials ([I]), (j2J) are determined from a given 5-matrix. In the 
previous version 3|, |J, [5[] we have restricted ourselves to the case of a two-channel system 
without a threshold. 

In this paper we attempt to extend the inverse scattering procedure to the two 

coupled channels with different thresholds: Ai = 0, A2 = A > 0. The channel wave numbers 
k a are related by 

k 2 a = k 2 -A a . (4) 

We measure the energy E in to, i. e. E = hue and e = q 2 /2, where q = pk. Thus, 
the elements of the Hamiltonian J\f x H (Af = 2N) matrix H in the basis ||ot) (f)^ , 
n = 0, . . . , N — 1, a = 1, 2} are written as 



Tjo.fi r J rp(ot) _|_ r P A \ _i_ \/a 13 



(5) 



Here, T^L = T%» m are the elements of the symmetric tridigonal matrix 

Tl n+1 = Ti +hn = -iyj(n + l) (n + e+i) 



(6) 



of the £th partial wave kinetic energy operator 

p 2 ( d 2 
~2 



dr 2 



(7) 



Recall that generally, with a finite order potential matrix it is possible to reproduce the 
scattering data only in finite energy interval. The S'-matrix (S(k)) features in a given interval 
k G [0, ko] determine the optimal combination of iV and p that is best suited to the task of 
the potential ([I]), (j2J) construction. 

Within the framework of the method the eigenvalues {Xj} and rows {Z^j}, {Ztfj} of the 
eigenvector matrix Z of the Hamiltonian matrix H are derived from the S'-matrix. As has 
been shown in Ref. [5j, the set {Xj, Zjyj, suffices to determine the Hamiltonian 

matrix H of the quasitridiagonal form 



/ a?> b? 
b? a? 
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(2) J2) ,(2) 



a 



b {2) a (Z > b 

u N-2 a N-2 U N-1 







(2) J2) 
N- 

a (2) / 

AT-l a N-l / 



(8) 



The present paper is organized as follows. In Sec. II we outline the multichannel J- 
matrix formalism. Sec. Ill is devoted to the generalization of the J-matrix inverse scattering 
formalism to the case of two coupled channels in the presence of the threshold. An example 
illustrated all steps of the presented inverse procedure is given in Sec. IV. Conclusions are 
drawn in Sec. V. The Appendix consists of a brief discussion of the J-matrix version of the 
two-channel Marchenko equations. 
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II. PRELIMINARIES 



The components ip( a/3 \k, r) of 2 x 2 matrix solution to the coupled radial Schrodinger 
equation, within the J-matrix formalism, are expanded in the oscillator function series [2] 

oo 

^)(fc,r) = ^^(x)cl^)(A;). (9) 



n=0 



For the potential ([T]), ([2]) the expansion coefficients c£^\k) are represented for n > N — 1 
as a combination 



Ho 




(10) 



of two linearly independent solutions 

Sn,M) = 



r(n+£+f) 



(11) 



Cn ,( g ) = /^liwa e v/2 F( _ n 

to the three-term recursion relation 



1/2, -£+1/2; q 1 



T n,n-idn-i(q) + T^ n d n (q) + Tl n+1 d n+1 (q) = —d n (q), n = 1, 2, 
which is the basis-set representation of the free Schrodinger equation. 



;i2) 



In turn, the component ipi 01 ^ of the bound eigenstate of energy E v = —^f^ expansion 



ft 2 *?. 



(<*/?), 



(13) 



n=0 



coefficients satisfy the boundary conditions 



d a/3) (is) = 5 aP ft\is), n>N-l, 



(14) 



The completeness relation of the scattering and bound states [7| can be transformed for 
the solutions 



Cn ^(^) Cn (&) 
Cn ' ( /c ) Cn ( & ) / 



(15) 



into 

oo 

2 



7T 



J dfcc n (fc)Pc m (fc) t + ^cMA 1/C M t = 15 n , T 







I being the 2x2 unit matrix. Here the matrices P and A are defined as 



where 



are the asymptotic normalization constants, i. e. the normalized bound-state 
function 

'^ 1) {t)^ 2 \t)\(m^ 
4 21 \r) 4 22) (r) j \M? ] 

asymptotically behaves as 



16] 



: 


^ 


1 




k Re(k 2 ) 




k 2 \k 2 \ J 





i Res S a ,(k) = ^J ^ + ^g + ^^q 



[ *8> > 


1 ~ ( 







The J-matrix expressions for the S'-matrix elements have the form 

£jV-l,£ 2 fe) ~~ ^22(e) ^jV-l.A^ife)) - ^12 (<0 Tn-X^N^N^nIM^N^M 2 ) 



D(+)(k) 
where 

D(+)(A;) = -Vn(e)T^ N C { +l( qi )) x 

x (c^fe) -^(6)^^(4.)) -n(e)Ti r i i 1 xv^i 1 (?i)4 H :ife)- 

The functions V a p(e) are defined by 

* e — a,- . * e — a,- r - : e — A, 

J= i 3=i J 3=1 J 

The eigenvalues Aj are thus the poles of V a p{e). 



III. DESCRIPTION OF THE METHOD 



Region y/A < k < ko 

In this energy region all channels are open, the S'-matrix is unitary. Thus the method 

0,9 

can be applied to calculate the eigenvalues Xj G [%-A, tt^o] an d corresponding eigenvector 
components Z N j, Z^-^. Firstly, the functions V a p{e) are defined by inverting (l2Tl) - (j23]) 
relative to V a p{e), 



fti(e) = m Me) = m VM @3(e) 



D(e) 



D(e) 



D{e) 



where 



(26) 



(27) 



C$ i (q 1 )-C£> i (q 1 )S ll (k)) [C£l 1/2 (q 2 ) - Cgl^fa) S 22 (k) 



»(+) 



(-) 



»(+) 



-1,JV 

ip 2 v^Si 2 (fc) 

T (l) T (2) 
1 N-1,N 1 N-1,N 



,(2) 



6,(e) 



(28) 



(29) 



D(e) = f C$ x (9i) - Cj&foi) 5 U (*)) (C^Cft) - 41(^2) S 22 (A;) 



»(+) 



»(-) 



»(+) 



CH{qi)C^l{q 2 )SUk). 



Then the eigenvalues A,- are the roots of 



D(e) = 0. 



(30) 



(31) 



In turn the residues of V a p(e) (126]) at the poles Xj determine Z^j, Zj^-j squared and the 
component products: 



Zjf d = ResVnie), 



(32) 



Zh, = ResV 22 (e), 



(33) 
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Z NJ Z MJ = ResV 12 (e). (34) 

e=\j 

Notice that the unitarity of the S-matrix provides that (for iV sufficiently large) the poles 
and residues of the functions V a p{^) (|26l) are real: 



Im{\,) = 0, Im UtesVapie)^ = 0, (35) 

and 

ResV aa {e) > 0. (36) 



Region k < y/A 

To define the functions V a p (1261) for this energies an analytic continuation of the S'- 
matrix to k 6 [0, v^A] would be carried out. To our knowledge there is no general method 
for doing this. However, the S'-matrix couldn't behave here in an arbitrary way. Generally, 
the use of non- unitary S-matrix can violate the conditions (|35|) . (|36|) . One way to meet the 
conditions is to restrict oneself to the using of the only open-channel (unitary) submatrix of 
the S-matrix. By this is meant that Eqs. (12S|) - (l3"0"|) reduce to 

^) = f# (37) 
D{e) 

e\(e) = (4i A (?i) " Sn{k)) /T^_ 1N , (38) 

D(e) = C { 'l ( gi ) - C$1 ( qi ) S n (k). (39) 

2 

Then the eigenvalues Xj < yA and the corresponding components Z N j of the eigenvectors 
are found from 

D(e) = (40) 

and 

Z? NJ = ResVi{e). (41) 

e=A J 

Finally, the components Zj^j are assumed to be zero. 
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Region k > fa 



o 



As noted before [5|, generally, a potential ([I]), ([2]) of finite rank can describe the scattering 
data only in a finite energy interval. In other words, the open-channel submatrix of the 
resulting S'-matrix coincides with the reference one on the interval [0, k ], in the region 
k > ho, however, such is not the case. On the other hand, it is necessary that there exist 
eigenvalues Xj > \k\. Since we have no prior knowledge of the S'-matrix for k > k , we 
couldn't draw on the Eqs. (I3~TT) - (I3"4"I) to evaluate the "external" parameters {Xj, Z N j, Z^-j}, 



\ . -> elk 2 

A 3 > 2 K Q . 



,5 



In the previous papers (see, e. g. [a]) the "external" parameters are obtained through 

the use of a standard fit to the data on the interval k G [0, ko\. In doing this the constraints 

AT AT AT 

E Z L = 1 ' ZXj = 1 ' E^^ = ' ( 42 ) 

3=1 3=1 3=1 

that follow from orthogonality of the eigenvector matrix Z, are allowed for the parameters 
Zmj, Ztfj. In the paper to this end we use the discrete version of the Marchenko equations 
js] generalized to the two coupled- channel case (see Appendix). Let us assume that, as in 
an example discussed below, only two largest eigenvalues Aj lie to the right of \k\. A set of 
six equations, that the "external" parameters {Xj, Znj, Zj^^}, j = N — 1, N are found from, 
must contain ( 142]) . In addition to this, equations can be derived using relations between the 

elements of the Hamiltonian matrix (jSJ) and its eigenvalues and eigenvectors. In particular, 

AT AT AT 

a N^l = E ^3 Z N,j, a N-l = E ^3 Z lf,ji U N~l = E Aj Z N j Zfifj- (43) 
3=1 3=1 j=l 

To obtain the quantities a^_ v ojv-i anc ^ u N-t the Marchenko equations (1A9I) . (1A10I) are 
solved for K^n , n = N — 2. Then Eqs. (TATll -l fATHll are used. 

Before we can use the Marchenko equations, we need to define the matrix S(k) in the 
integrands on the right-hand side of (1A6j) for < k < oo. For this purpose suppose that 
there exists k max > such that 

S(k) = I, k > k max . (44) 
Notice that the integrands in Eq. (IA6I) contain the solutions C n j(q) which are exponentially 



large as q — ► oo [Sj: 



C ne (q) ~ {-l) n+1 y/pn\ r(n + £ + 3/2)/vr q ~^+^) e i 2 /\ (45) 



q— too 



Thus the assumption (1441) provides the convergence of the integrals in (1A6j) . On the other 
hand, since the goal is to reproduce the scattering data in the interval k e [0, ho], it is 
natural to assume k max = k . Clearly, the open-channel submatrix of the matrix S(k) in the 
interval [0, ko] is determined by the scattering data. 

We also define the matrix S(k) in Eq. (1A6|) for energies where the second channel is 
closed. Notice that because of projection onto open channels (the matrix P in Eq. (1A6I) ) 
the contribution from the closed channel to the kernels Qnm reduces to the integrals over 
the interval [0, \/A] of the terms containing fn l \k) = — ^C^j 2 (q 2 ) -^J^ Si 2 (k) . In turn, for 
large n the solution C^j(i\q\) behaves as 9] 

C^}(i\q\) ~ %~ l Vp{n + t/2 + 3/4)"* e ^ q W n H+l. (46) 

( 21) 

Hence choosing A" large enough, the integrated terms involving /„ '(k) (n > N — 2) for 
k < \/~K can be made as small as we please. Thus for A" sufficiently large the approximation 

S 12 (k) = 0, k < VA (47) 

is acceptable for the matrix S(k) in Eq. (1A6j) . 

Let S^(k) denote the S'-matrix corresponding to the Hamiltonian matrix (jSJ) obtained by 
using the approximations (jHJ), fj47]) . Clearly, the properties (j4"4l) . f|4T|l of the matrix S(k) are 
not shared by the resulting matrix S^(fc). Thus, the following iteration procedure allow us 
to evaluate a contribution from the closed channel. At each iteration (i), in the integrands 
of the integrals over the interval [0, y/A] on the right side of ( 1A6I) that contain S\2(k) the 
element of the matrix S^ i_1 ^(A;) obtained at the previous step is used. The procedure is 
carried out until convergence is achieved. 



Bound states 

To every bound state with k u , A4u there corresponds the triplet {X u , Z Nt „, Zj^ jV }. 

The "bound" parameters {A^, Z N)U , Zj^ u } and |/c„, M?v \ M.f^ are related by 

DW(k„) = (48) 

and Eqs. ({TBI . Notice that \ u , Z N>V1 Z^^ v appear in the left members of Eqs. (T48|) and (|T8|) 
through the functions V a p fl25|) . In turn, k uj Ai^ are involved in the kernels Q^nm 
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(1A6I) of the Marchenko equations. Thus all of the unknown parameters {A-,-, Zjyj, Zj^-j}, 

2 

Xj > yA and {A„, Z N ^ V , Zj^^} are found by solving the system that consists of f|42l . fj43l) . 
(BHD and two out of four Eqs. (PSJ. 



IV. EXAMPLE 



As an example, we consider the S-matrix 



o _ {X — i k)(a 2 — b 2 + i a k — i a k 2 + k k 2 ) 
1 ' 'x + ik)(a 2 — b 2 — iak — iak 2 — kk 2 ) 



s = -lib ^kk 2 i^x 2 + A 2 - i k 2 ) (AQ) 
12 (x + ik)(a 2 — b 2 — iak — iak 2 — k k 2 ) ' 

g _ (V x 2 + A 2 — i k 2 )(a 2 — b 2 — iak + iak 2 + kk 2 ) 
(V x 2 + A 2 + i k 2 ) (a 2 — 6 2 — i ak — i ak 2 — k k 2 ) ' 



corresponding to a model s-wave 2x2 potential discussed in Ref. [10(] with a = —2, b = 0.6, 
x = 3. A threshold energy A = 10 is assumed in the second channel. (We take fi = fi = 1.) 
We found that the number = 5 of the basis functions ([3]) used in the expansion ([2]) and 
the oscillator radius quantity p = 0.495 are best suited to the task of a potential ([1]), (T5]) 
construction that reproduces the ^-matrix ()49l) in the interval k e [0, ko], &o = 6. 

The eigenvalues Xj, j = 4, . . . , 8 lie in the energy interval [yA, ^y 51 ] where the second 
channel is open. Then this eigenvalues and the corresponding eigenvector components are 
calculated from (!32l - (l34j) . The results are displayed in Table I. 

2 

The eigenvalues A-,- G [0, yA], j = 2, 3 (see Table I) and the corresponding components 
Zjvj, Zj^-j are found from the approximate Eqs. (140l - (|4TI) . Notice that if we use the ex- 
pression ( |49l) for the S'-matrix in ( l30l) . we obtain that the equation ( !3~T1) has not real roots 

A, < £a. 

„2 fc 2 

The calculations of the eigenvalues Xj > and the corresponding eigenvector com- 
ponents are carried out in combination with the description of the bound state that the 
system sustains at E b = —k 2 /2, k = 2.1946752413. The corresponding residues of the 
S'-matrix elements are 

Res S n (k) = -z26.7100700336, 

k ~ lK (50) 
ResS 12 (k) = zl8.1352046367. 

k=in 

The "bound" parameters {Ai, Z N ^, Z^i} and the "external" ones {Xj, Z N j, Zj^-j}, j = 
9, 10 are found by solving the system that consists of f|42l . fj43l . fj48|) . fl50|) . In the calculations 
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of <%_!, Wjv-i i n Eq. (Jl3"j) the approximations (jSj), (jUj) have been applied. The 

resulting parameters {Aj, Zjvj, ^Vj} are presented in Table I (a). Table 111(a) lists the 
corresponding Hamiltonian matrix (jSj). In Figs. 1-3 are shown the eigenphase shifts Si, S2 
and the mixing parameter e corresponding to the matrix S(k) in the integrands on the 
right-hand side of (1A6I) (thin solid line) and the resulting matrix S^(k) (dashed line). 

In order to evaluate the effect of the closed channel the iteration procedure described 
in the previous section is applied. The convergence of a^_ 1( o^-v u N-t 011 the left-hand 
side of (|43p depending on the number of iterations is displayed in Table II. The resulting 
parameters {Xj, Zjsrj, Zj^-j}, (j = 1, 9, 10) and the Hamiltonian matrix (jEJ) after 5 steps 
are listed in Tables 1(b) and Table 111(b) , respectively. The eigenphase shifts Si, S2 and 
the mixing parameter e corresponding to the matrix S^(fc) are plotted in Figs. 1-3 (dotted 
line) . 

A comparison of (a) and (b) results shows that, as we might expect, the off-diagonal 
part = V^ 21 ' of the interaction (pQ), ([2]) is most sensitive to the used approximations. 

The mixing parameter e behavior (see Figs. 3) reflects this dependence. Notice that the 
difference in the mixing parameters (in the interval [0, ho]) is within the accuracy of the 
method. 



V. CONCLUSION 



A generalization of the J-matrix inverse scattering approach to the case of two coupled 
channels with different threshold energies has been discussed. All the modification intro- 
duced in the inverse scheme 

aa 

(employed in the case of a two-channel system without 
threshold) appear to be plane and relatively obvious. An inverse procedure is proposed 
which focuss on reproducing the scattering data in a given energy interval. On the other 
hand, the procedure allows us to evaluate the contribution from the closed channel to the 
sought-for potential. 



Appendix A: THE J-MATRIX VERSION OF THE MARCHENKO EQUATIONS 

fn*^ (as well as C^) satisfies the three-term recursion relations ( Tl2l) . It follows herefrom, 
in view of the boundary conditions (fTUl) . (fl4j) and the Hamiltonian matrix quasitridiagonal 
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form (jSJ), that the coefficients ^^(k), Cn^\u) can be expressed as 

2N-n-2 



Ja/3) _ ST^ f 

n / j nm J i 



(7/3) 
m 



m=n 
7=1,2 



Besides, f lTUj) . fTi4"|) imply that K^fl N _ x = 5 a p. Further, from the expansion (IA1I) and 
completeness relations (fT6l) it follows that c ra is orthogonal to every f m for n < m, i. e. 



oo 



, m > n, 



fn(fc) 



where 

/^ 21) (fc) tf 2 \k) 

Inserting the expansion (1A1[) into (1A2j) and (fl6l) then yields 

2N-n-2 

/] K nn /Q„/ m = 0, m > 



rfV) 

f ( n\v) 



n. 



where 



Setting 



2N-n-2 



^ K nm Q mn K nn I, 



oo 

Qnm = ~ J dk f n (k) P f m (/c)^ + ^ ^ ^ A. u i^ , 



K r 



7^(11) if(12) 
~(21) _(22) 



we can rewrite (1A4j) in the form 

2V-n-2 



^ ^ Mnn'Qn'm Qnm; TYl ^> fl. 



n — n+1 



Inserting (|A8[) in (|A5[) gives 



K-77. 77, K- 



n n x Mi n 



2JV-n-2 



n -1 



m=n+l 
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Once M nm have been found by solving (1A9I) . K nn can be evaluated from (lAlOj) . Notice that 



in this case the Hamiltonian matrix quasitridiagonal form ([8]) implies that K nn = 0. 
A knowledge of K n a m finally furnishes the Hamiltonian matrix §8§ according to 



m m K (11) m K (12) K (21) 

A 1 ) rp(l) , f\in+lrp(l) _ -"-n, n -"- n -l, n rp 

«n — ±nn ~T ( n) J- n +l,n „(n)„(22) 1 i 



K {11) K {12) K (21) \ m 

-"-n-l. n J Si-l, n-l^n-l, n \ rp\i-j 



xl n-l,n— 1 xl n-l, n— 1 v n— 1, n-1 



(2) _ ^(2) I p A I V+ITW "'" A n „ + l rr (l) 

Un -intiT L\ t- (22 ) J n +l,n „(n) „(22) 



2 ~^7{2Z)~ n+1 , n „(ii) „(22) J-n+l.n" 



^(22) (12) (21) . 

xl n-l,n "■»> " "n-l, n | <JH ' 



„(22) „(11) „(22) ; x (l,n-l) 

^n-l.n-l n "n-l.n— 1 

^(22) 
-ill 



^(22) 

-"■71-1,71-1 



~(21) _(21) _(22) 



'(12) 7^(12) ~(11) 

■71,71 „(2) ^n-ljfi-l " 

7(22) _ -,(22) 7^(11) 



- £i n— 1, n— 1 xl n— 1, n— 1 Ji n— 1, n— 1 



(All) 



1 (1) 
^™ = 7^(11) ^7l,n-l> (^12) 

71—1, 71—1 



(A13) 



(2) /v n , n ( 2 ) 

7i „f99^ 7i, n— 1> y-^-^^J 



_ -^n,«+l T (l) J> -n-l,n ( 2 ) 

"ti ( n) ^n+l,n (ii) (22) J 7i,n-1) V-^ LO ) 

J^n,n J^n,n •"■n-l.rv-l 



7>W ~-n-i, Ti-i —, (X ) (AlK, 

n,n-l „f29.1 „/11l J ti, Ti-i' l^- ±u J 
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Table I: The eigenvalues and the eigenvector components of the Hamiltonian matrix ([8]). 



3 




A, 


Z N ,j 




i 
i 


(a) 


-0.58970264603 


-0.012333185006 


0.00018915595025 


(b) 


-0.58970263847 


-0.012333362733 


0.00019663447693 


2 




0.40533438179 


0.18715853083 





3 




0.78492505414 


0.090561490976 





4 




1.5244162751 


0.32069508364 


0.068124022185 


5 




1.6473065387 


0.14429242884 


-0.17954224586 


6 




2.6728511248 


0.056915520013 


0.34337676573 


7 




3.2807264076 


0.50263953170 


-0.052996407816 


8 




4.3839439596 


0.050257472502 


0.49798274116 


9 


(a) 


6.3135416565 


0.75676372010 


-0.062001763734 


(b) 


6.3156804888 


0.75741514694 


-0.048160786961 


10 


(a) 


7.5694793828 


-0.043006216886 


-0.76848969549 




(b) 


7.5651078745 


-0.029379475185 


-0.76948110211 
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Table II: Convergence of a N _i, <%_i> u N-l on the left-hand side of (fl3j) . i is the number of 
iterations. 

(1) (2) 

I a N-\ a N-l U N _i 

4.689928491 5.966326902 0.0191266184 

1 4.689911469 5.965705556 0.0071475853 

2 4.689912701 5.965663226 0.0059643414 

3 4.689912839 5.965658934 0.0058474387 

4 4.689912852 5.965658510 0.0058359055 

5 4.689912854 5.965658464 0.0058347978 
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Table III: The elements of the Hamiltonian matrix (JSj) . 



n 


a (1) 


<Jn 


(2) 

a> n 


,(2) 
"n 


Un 




(a) 





0.1775666649 




0.06753394526 




0.6517927001 




1 


1.212316541 


-0.3078456533 


2.378081035 - 


0.01876213797 


0.03459618331 


0.1112824167 


2 


2.233751540 


-0.8637975100 


3.353299231 


-0.7967451984 


0.02226447163 


-0.06006980906 


3 


3.357567064 


-1.370400686 


4.556450725 


-1.319787140 


0.04861624566 


-0.07835689367 


4 


4.689928491 


-2.007758312 


5.966326902 


-2.032477517 


0.0191266184 


-0.06589047115 


(b) 





0.2553490884 




-0.01023393318 




0.6404835124 




1 


1.212396865 


-0.3006058502 


2.378394045 - 


■0.01965488685 


0.03937726990 


0.1295217295 


2 


2.233837107 


-0.8638658833 


3.353460226 


-0.7967777718 


0.02976290436 


-0.06151639336 


3 


3.357902488 


-1.370738931 


4.553912267 


-1.319081324 


0.04174355881 


-0.07598084688 


4 


4.689912854 


-2.008020573 


5.965658464 


-2.031439617 


0.0058347978 


-0.04479790470 
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Figure 1: The eigenphase shift 5\ corresponding to the matrix S(k) on the right-hand side of (1A6 
(thin solid line) and the resulting matrices: S^(k) (dashed line) and S^ 5 \k) (dotted line). 
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Figure 2: The eigenphase shift 62- See Fig. 1 for details. 
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Figure 3: The mixing parameter e. See Fig. 1 for details. 
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